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We report on the extraction of the structure functions F?, and AxF^ — xFg — iF 3 " from CCFR i/ M -Fe 
and /7 M -Fe differential cross sections. The extraction is performed in a physics model independent 
(PMI) way. This first measurement of AxF$, which is useful in testing models of heavy charm 
production, is higher than current theoretical predictions. The ratio of the Fi (PMI) values measured 
in z/ M and fi scattering is in agreement (within 5%) with the predictions of Next-to-Leading- Order 
parton distribution functions (NLO PDFS) using massive charm production schemes, thus resolving 
the long-standing discrepancy between the two sets of data. 
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Deep inelastic lepton-nucleon scattering experiments 
have been used to determine the quark distributions in 
the nucleoli. However, the quark distributions deter- 
mined from muon H] and neutrino [0] experiments were 
found to be different at small values of Bjorken x, be- 
cause of a disagreement in the extracted structure func- 
tions. In this Letter, we report on a measurement of 
differential cross sections and structure functions from 
CCFR z^-Fe and v^-Fe data. The neutrino-muon differ- 
ence is resolved by extracting the v^ structure functions 
in a physics model independent (PMI) way. We also re- 
port on the first measurement of AxFj, = xF% — xF£, 
which is used to test models of heavy charm production. 

The sum of v^ and F M differential cross sections for 
charged current interactions on an isoscalar target is re- 
lated to the structure functions as follows: 
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Here Gf is the Fermi weak coupling constant, M is 
the nucleon mass, E v is the incident neutrino energy, 
the scaling variable y — EhjE v is the fractional energy 
transferred to the hadronic vertex, Eh is the final state 
hadronic energy, and e ~ 2(1 — y)/(l + (1 — y) 2 ) is the 



polarization of the virtual W boson. The structure func- 
tion 2xF\ is expressed in terms of F-i by 2xF\(x, Q 2 ) = 
F2(x,Q 2 ) x i+r(xQ2) ; where Q 2 is the square of the 
four-momentum transfer to the nucleon, x = Q 2 /2MEh 
(the Bjorken scaling variable) is the fractional momen- 
tum carried by the struck quark, and R — — is the ratio 
of the cross-sections of longitudinally- to transversely- 
polarized W bosons. 

A similar equation for the case of muon scattering re- 
lates the cross sections to the structure functions. How- 
ever, there are significant differences originating from the 
scattering on strange (s) and charm (c) quarks. The 
AxF% term, which in leading order ~ Ax(s — c), is not 
present in the fi scattering case. In addition, in a charged 
current v^ interaction involving s (or c) quarks, there is a 
threshold suppression originating from the production of 
heavy c quarks in the final state. For /i scattering, while 
there is no suppression for scattering from s quarks, there 
is larger suppression when scattering from c quarks since 
there are two heavy quarks (c and c) in the final state. 

In previous analyses []2j of v^ L data, light flavor univer- 
sal physics model dependent (PMD) structure functions 
were extracted by applying a slow rescaling correction to 
correct for the charm mass suppression in the final state. 



In addition, the Acci 7 ^ term (used as input in the extrac- 
tion) was calculated from a leading order charm produc- 
tion model. Recent calculations ||-||] indicate that there 
are large theoretical uncertainties in the charm produc- 
tion modeling for both AxF^, and the slow rescaling cor- 
rections. Therefore, in the new analysis reported here, 
slow rescaling corrections are not applied, and Aa;^ and 
F2 are extracted from two-parameter fits to the data. We 
compare the values of AxF^ to various charm produc- 
tion models. The extracted physics model independent 
(PMI) values for F% are then compared with F% within 
the framework of NLO models for charm production. 

The CCFR experiment collected data using the Fermi- 
lab Tevatron Quad- Triplet wide-band v^ and v ^ beam. 
The CCFR detector || consists of a steel-scintillator 
target calorimeter instrumented with drift chambers, 
followed by a toroidally magnetized muon spectrom- 
eter. The hadron energy resolution is AEh/Eh = 
0.85/v / -E-h(GeV), and the muon momentum resolution is 
Ap^/pp, — 0.11. By measuring the hadronic energy (Eh), 
muon momentum (py), and muon angle (0^), we con- 
struct three independent kinematic variables x, Q 2 , and 
y. The relative flux at different energies, obtained from 
the events with low hadron energy (Eh < 20 Ge V), is nor- 
malized so that the neutrino total cross section equals the 
world_average a uN /E = (0.677±0.014) x 10" 38 cm 2 /GeV 
and a VN /a uN = 0.499±0.005 @. After fiducial and kine- 
matic cuts (p^ > 15 GeV, 0,t < 0.150, E h > 10 GeV, and 
30 GeV < E v < 360 GeV), the data sample consists of 
1,030,000 v^ and 179,000 v^ events. Dimuon events are 
removed because of the ambiguous identification of the 
leading muon for high-y events. 

The raw differential cross sections per nucleon on iron 
are determined in bins of x, y, and E v (0.01 < x < 0.65, 
0.05 < y < 0.95, and 30 < E„ < 360 GeV). Figure 
[j] shows typical differential cross sections at E v = 150 
GeV (complete tables are available J7J). For all en- 
ergies, the cross sections are in good agreement with 
NLO PDFs (with massive charm production schemes 
e.g., MRST99 @ or CTEQ4HQ §). The dashed lines 
shows the predictions from the Thome and Roberts Vari- 
able Flavor Scheme (TR-VFS) @ QCD calculation using 
MRST99 extended Q PDFs. This calculation includes 
an improved treatment of massive charm production. 
The QCD predictions, which are on free neutrons and 
protons, are corrected for nuclear 0] , higher twist pl|^2| , 
and radiative effects jl3| . Also shown are the predictions 
from a CCFR leading order (LO) QCD inspired fit used 
for calculation of acceptance and resolution smearing cor- 
rections. As expected from QCD, the CCFR cross section 
data exhibit a quadratic y dependence at small x for v^ 
and v ' n, and a flat y distribution at high x for v^. 

Next, the raw cross sections are corrected for elec- 
troweak radiative effects [O, the W boson propagator, 
and for the 5.67% non-isoscalar excess of neutrons over 
protons in iron (only important at high x). Values of 



do-VEdxdy : 150 GeV 




7 [ v ] y [ v ] 

FIG. 1. Typical raw differential cross sections at E v = 150 
GeV (both statistical and systematic errors are included). 
The CCFR data are in good agreement with the NLO 
TR-VFS QCD calculation using MRST99 PDFs (dashed 
line). The solid line is a CCFR LO QCD inspired fit. 



ASF3 and F% are extracted from the sums of the cor- 
rected v^-Fe and F M -Fe differential cross sections accord- 
ing to Eq. (1). However, it is challenging to fit AxF§, 
R, and 2a;Fi using the y distribution at a given x and 
Q 2 because of the strong correlation between the AxF^ 
and R terms, unless the full range of y is covered by the 
data. Covering this range (especially the high y region) 
is difficult because of the low acceptance. Therefore, we 
restrict the analysis to two-parameter fits. 

Our strategy is to fit AxF% and F^ (or equivalently 
2a; Fi) for x < 0.1 where the AxF-j contribution is rela- 
tively large, while constraining R using the R^otM 01 
QCD-inspired empirical fit to all available electron and 



muon scattering data. The R^rid n ^ ^ s a ^ so m good 
agreement with NMC i? M data JJJ at low x, and with the 
most recent theoretical prediction |12| for R (a NNLO 
QCD calculation including target mass effects). For 
x < 0.1, R in neutrino scattering is expected to be some- 
what larger than R for muon scattering because of the 
production of massive charm quarks in the final state. A 
correction for this difference is applied to R^ r i d using a 
leading order slow rescaling model to obtain an effective 
R for neutrino scattering, R^tr- The difference between 

R-worid an d Reff ^ s use< ^ as & systematic error. Because 
of the positive correlation between R and Aa^i 7 ^, the ex- 
tracted values of i*2 are rather insensitive to the input 
R. In contrast, the extracted values of AxF% are sensi- 
tive to the assumed value of R, which is reflected in a 
larger systematic error. The values of Aa;F3 are sensi- 
tive to the energy dependence of the neutrino flux (~ y 
dependence), but are insensitive to the absolute normal- 
ization. The uncertainty on the flux shape is estimated 
by using the constraint that F^ and a;F3 should be flat 
over y (or E v ) for each x and Q 2 bin. 

Because of the limited statistics, we use large bins in 
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FIG. 2. AxF$ data as a function of x compared 

with various schemes for massive charm production:(left) 
TR-VFS (MRST99), ACOT-VFS (CTEQ4HQ), FFS(GRV94), 
and the CCFR-LO (a leading order model with a slow rescal- 
ing correction) ; (right) Sensitivity of the theoretical calcula- 
tions to the choice of scale. 



Q 2 in the extraction of AxF^ with bin centering correc- 
tions applied using the NLO TR-VFS calculation Q] with 
the MRST99 PDFs. Figure | (left) shows the extracted 
values of AxF^ as a function of x, including both sta- 
tistical and systematic errors, compared to various the- 
oretical methods for modeling heavy charm production 
within a QCD framework. The three-flavor Fixed Fla- 
vor Scheme (FFS) [[15[ assumes that there is no intrin- 
sic charm in the nucleon, and that all scattering from c 
quarks occurs via the gluon-fusion diagram. The concept 
behind the Variable Flavor Scheme (VFS) proposed by 
ACOT H[16| is that at low scale, li, one uses the three- 
flavor FFS scheme, while above some scale, an intrinsic 
charm sea (which is evolved from zero) is introduced. 
The concept of the TR-VFS scheme 0] is similar, except 
that at intermediate scale it interpolates smoothly be- 
tween the two regions. Both the FFS and VFS schemes 
have been implemented by KLS B, ACOT, and Kret- 
zcr JlTj ] - The last two implementations agree with each 
other, but not with KLS (there was a mistake in the KLS 
calculation) fl8|| . 

Shown are the predictions from the TR-VFS scheme 
(implemented with MRST99 PDFs and the suggested 
scale li = Q) and along with the predictions from two 
other other NLO calculations, ACOT-VFS (implemented 
with CTEQ4HQ PDFs and their recent suggested scale 
li = m c for Q < m c , and li 2 = m c 2 + 0.5Q 2 (l — m c 2 /Q 2 ) n 
for Q < m c with n = 2 ||), and the FFS (implemented 
with the GRV94 |l^] PDFs and their recommended scale 
li = 2m c ). Also shown are the predictions from the lead- 
ing order QCD fit to the CCFR dimuon H data. 
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FIG. 3. The ratio of F 2 " (PMI) data divided by F 2 M 
(NMC or BCDMS) or F| (SLAC). Both statistical and sys- 
tematic errors are included. Also shown are the predictions of 
the TR-VFS (MRST99), ACOT-VFS (CTEQ4HQ) and FFS 
(GRV94) heavy flavor calculations. 

Figure || (right) shows the sensitivity to the choice of 
scale. For example, the data do not favor the choice of 
scale, li = 2Pt max in the ACOT-VFS calculation with 
CTEQ4HQ PDFs. This high scale (originally suggested 
by ACOT and used in the CCFR dimuon analysis [§§) 
implies that the calculation is in the four-flavor region 
even at x — 0.015 and Q 2 = 1.0 GeV 2 (and yields large 
negative result). With reasonable choices of scale, all the 
theoretical models yield similar results. However, at low 
Q 2 , our AxFs data is higher than all of the theoretical 
models. The difference between data and theory may be 
due to an underestimate of the strange sea at low Q 2 , or 
from missing NNLO terms. The question of the strange 
sea would be addresed by a global NLO analysis which 
combines the neutrino data for dimuons, AxF^ and F, 
with Fg" e . 

As discussed above, values of F 2 for x < 0.1 are ex- 
tracted from two-parameter fits to the y distributions. 
In the x > 0.1 region, the contribution from AxF^ 
is small and the extracted values of F 2 are insensitive 
to AxF^. Therefore, we extract values of F 2 with an 
input value of R and with AxF^ constrained to the 
TR-VFS (MRST99) predictions. Fig. | shows the ra- 
tio of our F% (PMI) measurements divided by (18/5)F 2 M 
(NMC § or BCDMS @) or (18/5)F| (SLAC @) mea- 
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surements |23 . The overall normalization errors of 2% 
(CCFR), and 2.5% (NMC) are not shown. Within 5%, 
the ratio is in agreement with the predictions of the 



TR-VFS (MRST99), ACOT-VFS (CTEQ4HQ), andFFS 
(GRV94) calculations [||. 

In the calculation of the theoretical predictions, we 
have also included corrections for nuclear effects S. As 
mentioned earlier, the extracted values of Fq, from the 
two-parameter fits are insensitive to R. For example, if 
we perform simultaneous two-parameter fits to Fi and 
R (while keeping AxF 3 fixed to the TR-VFS (MRST99) 
values), the extracted R values at x = 0.01 are smaller 
than R"ff, but F 2 changes by only 2 ~ 3%. 

In the previous analysis of the CCFR data, the 
ratio of extracted values of F% (PMD) data divided by 
(18/5)F 2 M (NMC) at the lowest x = 0.015 and Q 2 bin 
were 20% higher than the predictions of the light-flavor 
PDFs such as MRSR2 || or CTEQ4M § (see Fig. | ). 
About 10% of the difference originates from having used a 
leading order model for AxF% versus using our new mea- 
surement. Another 6% originates from having used the 
leading order slow rescaling corrections, instead of using 
NLO massive charm production models. The remain- 
ing 3% originates from improved modeling of the low Q 2 
PDFs (which changes the radiative corrections and the 
overall absolute normalization to the total neutrino cross 
sections). For higher Q 2 at x = 0.015, and for the next 
two x bins (x = 0.045 and 0.08), the smaller difference 
between the PMI and PMD results is due to equal con- 
tributions from the AxF% and the difference in the slow 
rescaling corrections. For the higher x bins (x > 0.1), 
the contribution of A2F3 is small, and the slow rescaling 
corrections in the leading order model are the same as 
those with the NLO theories. Therefore, the NMC and 
CCFR data are in agreement at large x whether PMI or 
PMD structure functions are used in the comparison. 

In conclusion, the ratio of F% (PMI) values measured 
in neutrino-iron and muon-deuterium scattering are in 
agreement with the predictions of Next-to-Leading-Order 
PDFs (using massive charm production schemes), thus 
resolving the long-standing discrepancy between the two 
sets of data. The first measurement of AxF^ is higher 
than current theoretical predictions. 
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